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ABSTRACT 

We study the evolution of the luminosity-to-halo mass relation of Luminous Red Galaxies (LRGs). We select a sample 
of 52 000 LOWZ and GMASS LRGs from the Baryon Oscillation Spectroscopic Survey (BOSS) SDSS-DRIO in the ~450 
deg^ that overlaps with imaging data from the second Red-sequence Gluster Survey (RGS2), group them into bins of 
absolute magnitude and redshift and measure their weak leasing signals. The source redshift distribution has a median 
of 0.7, which allows us to study the leasing signal as a function of lens redshift. We interpret the leasing signal using 
a halo model, from which we obtain the halo masses as well as the normalisations of the mass-concentration relations. 
We find that the concentration of haloes that host LRGs is consistent with dark matter only simulations once we allow 
for miscentering or satellites in the modelling. The slope of the luminosity-to-halo mass relation has a typical value of 
1.4 and does not change with redshift, but we do find evidence for a change in amplitude: the average halo mass of 
LOWZ galaxies increases by 251^®% between z = 0.36 and 0.22 to an average value of (6.43T0.52) x 10 If we 

extend the redshift range using the GMASS galaxies and assume that they are the progenitors of the LOWZ sample, 
we find that the average mass of LRGs increases by 80ji5o% between z = 0.6 and 0.2. 


Key words, gravitational lensing - dark matter haloes 


1. Introduction 

Hierarchical models of structure formation predict that 
galaxies form in small dark matter haloes, which sub- 
sequently cluinp tog ether and merge into larger ones 
(|White fc ReesI 119781) . At large scales, the evolution of 
structure is mainly determined by the properties of dark 
matter and dark energy. However, at smaller, galactic 
scales, baryonic physics cannot be ignored. Processes such 
as supernova and AGN feedback impact the relation be¬ 
tween the observable (baryonic) properties of galaxies and 
their dark matter haloes. Hence by measuring these re¬ 
lations, we gain insight into the processes that affected 
them. Studying this with numerical simulations is noto¬ 
riously difficult, although in recent years this field has 
rapidly a dvanced thro ugh the use of semi-analytic mod¬ 
els ( e.g. iBau ^ 120061) and hydrodynamical sim ulations 
(e.g. IVogelsberger et al]l20l4 ISchave et al.ll2015l) . To test 
these simulations and guide them with further input, we 
need observations of the relation between the properties 
of galaxies and their dark matter haloes. This is also cru¬ 
cial for understanding the effect of baryonic physics on the 
dark matter power sp ectrum (e.g Ivan Daalen et all 120111: 
ISemboloni et al.l 120111) , which is the main observable in 
weak lensing studies that aim to extract cosm ological pa¬ 
rameters, such as Euclid (|Laureiis et al.l[20Tl^ . 

The properties of dark matter haloes around galaxies 
can be studied with weak gravitational lensing. As the pho¬ 
tons emitted by distant galaxies traverse the Universe, they 
are deflected due to the curvature of space around interven¬ 
ing mass inhomogeneities in the foreground. Consequently, 
the observed shapes of these background galaxies slightly 


deform, a distortion that can be reliably measured out to 
proj ected separations of tens of Mpcs around the lenses 
(e.g. lMandelbaum et al.l[2?)T^ . Since this completely covers 
the regime where the dark matter halo of any lens domi¬ 
nates, weak gravitational lensing offers an excellent tool to 
measure halo masses. The weak lensing signal of individual 
galaxies is too noisy to be detected, but by averaging the 
signal of many lenses of similar observable properties, e.g. 
in a certain luminosity range, we can learn about the aver¬ 
age halo properties of such lens samples. 


The relation between the properties of galaxies and their 
dark matter haloes has been studied before with weak lens- 


ing (e.g. Hoekstra et al.l 2005 

: iMandelbaum et akl l2006b 

Li et al.ll2009l: van Uitert et al. 

I 2 OII: Brimioulle et al.ll2013 

Velander et al.ll2014D. but most of these studies focused 


on lenses at a limited redshift range. However, to study 
how galaxies evolve, one would like to measure how the 
luminosity-to-halo mass relation depends on lookback time. 
Recent imaging surveys such as the Canada-France-Hawaii 
Telescope Survey (CFHTLS) and the second Red-sequence 
Cluster Survey (RCS2) contain sufficient statistical power 
to enable such studies. Redshift dependent constraints that 
are derived in a homogeneous way, as is done in this study, 
are particularly useful for numerical simulations, as they 
can potentially disentangle degeneracies among the model 
parameters and limit the space for fine-tuning to match 
low-redshift obs ervations (for an example, see Figure 23 in 
iGuo et ^120111) . 

In this work, we study a particular type of galaxies: 
Luminous Red Galaxies (LRGs). They form an interest¬ 
ing subsample of the total population of galaxies, as they 
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trace the highest density peaks in the Universe. These 
galaxies are thought to have formed around z^2 during 
a relatively short and intense period of star formation, af¬ 
ter which the formation of stars practically halted. Their 
luminosity evolution can therefore be approximately de¬ 
scribed by ‘passive evolution’, the evoluti on of a stellar pop- 
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large the average impact is on the luminosity evolution, 
compared to the pure passive evolution scenario, is not 
clear. However, for massive and luminous LRGs, the lu¬ 
minosity evolution is thought to be well understood. 

Also from an observational perspective, LRGs are ad¬ 
vantageous to study. They are easily selected in multi-band 
optical data-sets and their redshifts can be relatively easil y 
determined using the 4000A break (|Eisenstein et al.ll2001f) . 
More than a million LRGs have been observed spectroscopi¬ 
cally as part of the Baryon O scillation Spectroscopic Survey 
fBOSS: Ibawson et al.ll201^ . forming the LOWZ sample, 
which targets z < 0.4 galaxies, and the CMASS sample, 
which targets 0.4 < z < 0.7 galaxies. From a weak leas¬ 
ing perspective, the advantage of LRGs is that they are 
massive and therefore produce a large leasing signal that 
can be measured up to relatively high redshift. The overlap 
between the BOSS survey and the RCS2 therefore offers 
a perfect combined dataset to study the evolution of the 
luminosity-to-halo mass relation of LRGs. 

The outline of this work is as follows. In Section [5] we 
describe the data that we use in this work, how we com¬ 
pute the luminosities and how we perform the leasing anal¬ 
ysis. We interpret the leasing measurements with the halo 
model, which we describe in Section [S] The evolution of the 
luminosity-to-halo mass relation is presented and discussed 
in Section 01 The mass-concentration relation is discussed 
in Section [5l We conclude in Section [6] Un less stated oth¬ 
erwis e, we assume a WMAP7 cosmology ([Komatsu et al.l 
IMlh with 0-8 = 0.8, Ha = 0.73, Hm = 0.27, Hb = 0.046, 
and /i7o = Ho/70 km s“^ Mpc“^ with Hq the Hubble con¬ 
stant; all distances are quoted in physical (rather than co¬ 
moving) units; and all apparent magnitudes have been cor¬ 
rected using the dust maps from ISchleerel et all 


2. Data analysis 

In this work w e make use of d ata from the tenth data 
release (DRIO; lAhn et ^ 2014) fro m the Sloan Digital 
Sky Survey ISDSS: lYork et ah jpOOf) and from the secon d 


Red-s equence Cluster Survey (RCS 2; Gilbank_e£jb 201T). 
As in Ivan Uitert et all (l201lL[20T^ : ICacciato et aLr ( 2014 ). 


we take advantage that the SDSS contains more ancillary 
data on galaxies than is available in the RCS2 due to its 
photometry in five optical band and its spectroscopy for 
over a million of galaxies. However, the RCS2 imaging is 
^2 magnitudes deeper and achieved a median seeing of 


approximately 0.7", compared to 1.2" for SDSS, making 
the RCS2 better suited for a weak lensing analysis of lenses 
at higher redshifts. The total overlap between the RCS2 
and the DRIO amounts to roughly 450 square degrees. 


A first combined analysis of the overla 

p between the 

ninth data release of SDSS IDR.O: |Ahn et al. 

|2012|) and the 

R.CS2 was presented in ICacciato et al.l (201 

4), where the 


lensing signal of the DR9 galaxies with spectroscopy was 
studied using RCS2 galaxies as sources. In that work we 
did not study the redshift evolution of the lensing signal, 
as, in contrary to the current work, we studied a mixed 
sample of early- and late-type galaxies, whose combined 
luminosity evolution was not well understood. Also, in 
DR9, the number of (high-redshift) BOSS spectra was 
roughly half of that in DRIO. 


2.1. Lens sample 

We use a subset of the total sample of overlapping DRIO 
galaxies with spectroscopy as our lenses, i.e. only the LRGs. 
We select all galaxies that have been targeted as part of 
BOSS. These are selected from the SDSS catalogues by 
requiring Q 

• BOSS_TARGETl && 2° 

• SPECPRIMARY == 1 

• ZWARNING.NOQSO == 0 

• TILEID >= 10324 

for the LOWZ sample, and 

• BOSS_TARGETl && 2^ 

• SPECPRIMARY == 1 

• ZWARNING.NOQSO == 0 

• (CHUNK != “bossl”) && (CHUNK != “boss2”) 

• *fib2 < 21.5 

for the CMASS (high-z) sample. Additionally, we select 
all objects with reliable spectroscopy from the SDSS cata¬ 
logues that satisfy the BOSS LOWZ target selection cuts: 

• \ir-i)-ig- r)/^ - 0.18| < 0.2 

• r < 13.5 -f [0.7 X (g - r) -f 1.2 x ((r - i) - 0.18)]/3 

• 16 < r < 19.6 

• SCIENCEPRIMARY==1 

• ZWARNING.NOQSO == 0 

• ^spec ^ 0.01 

where g, r and i indicate model magnitudes and r cmodel 
magnitudes. Note that we replaced the BOSS selection cri¬ 
terion Tpsf — Tcmod > 0.3 with Zspec > 0.01 to ensure that 
we have no stars. Finally, we also selected all objects that 
satisfied the CMASS selection cuts, but we found that all 
objects were already targeted and labeled as being BOSS 
galaxies, and it therefore did not increase the lens sample. 

Even though the LOWZ and CMASS samples mainly 
consist of LRGs, the populations diffe r due to the d i fferent 
colour and magnitude selection cuts. iToieiro et al.l (I2012D 
study which fraction of the CMASS LRGs are progenitors 
of the LOWZ sample, and find that this strongly depends 
on absolute magnitude, with the highest fractions found for 
the most luminous objects. A second but weaker trend is 

^ http://www.sdss3.org/dr9/algorithms/boss_galaxy_ts .php 
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the luminosities of our LRGs, corrected for the redshift of 
their spectra through the passbands (i.e. the fc-correction). 
We compute the fc-co r rection using the KCORREC T v4_2 code 
(|Blanton et al.l 1200,11 : iBlanton fc Roweii l2007li . where we 
use the tt, g, r, i and 2 model magnitudes and the spec¬ 
troscopic redshift as input. Furthermore, we correct for the 
intrinsic evolution of the luminosities (the e-correction), ac¬ 
counting for the difference between the observer-frame ab¬ 
solute magnitude of a galaxy with and without an evolving 
spectrum. 

The luminosities of LRGs are thought to evolve pas¬ 
sively, which can be modeled using a stellar population syn¬ 
thesis code. We make use of one of the p ublicly available 
codes, GALAXEV (iBruzual fc Charl^ 2003fl. i n the default 


configuration, i.e adopting a lGhabrier ( 200.111 IMF and us- 


Fig. 1. The SDSS g — r and r — i colours versus redshift 
for the LRG sample used in this work. The solid cyan 
li ne indicates the me d ian. T he other lines show a range 
of IBruzual &: Gharlotl l|2003ll SSP models, for various for¬ 
mation times (different colours) and metallicities (different 
line-styles) as indicated in the top-left of the figure. 


found with rest-frame colour. Therefore, we choose to anal¬ 
yse the LOWZ and the CMASS samples separately. We do 
investigate what we can conclude about the evolution of the 


ing the Padova 1994 tracks for the stellar evolution. We 
compute a range of instantaneous-burst models, where we 
vary the formation time and the metallicity. In Figure [U 
we show the evolution of the g — r and r — i colours of these 
models, together with the observed colours of the LRGs. 
The set of models that describe the data best are those 
that assume a metallicity oi Z = 0.02 (Zq). However, at 
2; < 0.4, the observed g — r colours are a bit too red, and 
at 0.4 < z < 0 7 the r — i colours are somewhat too red. 
iMaraston et al.l (1200911 improved the modelling by includ¬ 
ing a very low metallicity component to t he mod e l that 
consisted of 3% in mass, and by using the iPicklesI (|l998ll 
empirical spectral library instead of the theoretical one. 
However, below a redshift of ~0.5 the evolution correc¬ 
tion is fairly insensitive to the details of the modelling (see 
Figure [2]), while at higher redshi f ts it i s not clear whether 
the changes from iMaraston et al.l (|2009ll improve the match 
due to the low number of objects at this redshift range used 
in that work. As we will discuss below, our results do not 
critically depend on the choice of the model, hence we do 
not deem it ne c essary to include the improvements from 
IMaraston et al.l (l2009ll . 

In Figure[2]we show the fc-correction that these GALAXEV 
models predict, together with the fc-correction for the LRGs 
that have been computed using KCORRECT. We find that at 
2; < 0.4, the Z = 0.02 tracks agree well, but at higher red- 
shifts the fc-correction values from KCORRECT are somewhat 
lower than the GALAXEV model. In fact, the agreement at 
0.4 < z < 0.7 with the Z = 0.008 metallicity models is 
remarkable, but the validity of these models for our LRGs 
at low redshift is excluded based on the colour evolution in 
Figured] However, at z > 0.4 the LRGs show an increasing 
scatter in their colours and become more compatible with 
the Z = 0.008 models. 

We show the luminosity evolution of some of the 
GALAXEV models in Figure[2| We only show the models with 
Z = 0.02 and Z = 0.008, as the models with different metal¬ 
licities are excluded based on their colour evolution and k- 
corrections. Also, we only show models with a formation 
time of 10, 11 and 12 Gyrs, as most previous works on the 


luminosity-to-halo mass relation of LRGs if we consider the 

time 

in this range fe.g 

Wake et alJ 20061: 

Maraston et al. 

GMASS sample as progenitors of the LOWZ LRGs. 

2009 

:lBanerii et al.ll201C 

: iCarson & Nicholl 

201Ct iLiu et al. 


2012fl . For our nominal luminosity evolution correction, we 


2.1.1. Luminosities 

In order to study how the average halo mass of LRGs 
evolves, we want to compare LRGs at low redshifts to their 
predecessors at higher redshifts. Hence we need to obtain 


adopt the Z = 0.02 model that formed 11 Gyrs ago (at red¬ 
shift 0). As none of the models exactly captures the trends 
in Figured) and [21 the evolution correction we use may have 
a small bias. However, we have also tried different evolution 
corrections, with corresponding models that broadly cover 
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Redshift 

Fig. 2. {top panel:) The fc-correction as a function of red- 
shift. The gray scale shows the /c-corrections from the 
KCORRECT code, the solid cyan line indicates the median 
and the other colo ured lines show the fc-corr ection as pre¬ 
dicted by a range of iBruzual fc Charlotl (l2003ll SSP models, 
for various formation times (different colours) and metal- 
licities (different line-styles) as indicated in the top-left of 
the figure, {bottom panel:) The luminosity evolution cor¬ 
rection as a function of redshift. For c larity, we only show 
a few of the iBruzual fc Charlotl (|2003ll model predictions. 
The thick green dot-dot-dashed line shows the correction 
we have used in this work, which is based on the Z = 0.02 
instantaneous-burst model that formed llGyrs ago. 


the observed colour evolution and fc-correction values. We 
detail on this test in Section [dTj the main result is that our 
results do not change significantly. This suggests that the 
systematic bias in the luminosities caused by an incorrect 
evolution correction is likely insignificant for this work. 


LRGs have formed over a certain range of time and with 
some range of metallicities. Hence their actual luminosity 
evolution corrections may have some scatter compared to 
our nominal correction, as we found that the luminosity 
evolution correction is increasingly sensitive with redshift 
to these parameters. If this scatter is random with respect 
to our nominal correction, this causes an Eddington bias, as 
lenses are preferentially scattered to where there are fewer 
of them. In Appendix [21 we estimate the impact this may 
have on our masses estimates. We hnd that it is signih- 
cantly smaller than our statistical errors and can be safely 
ignored. 

In Figure |3l we show the distribution of absolute mag¬ 
nitudes after including the fc-correction and the (fc + e)- 
correction. In the range 0.15 < z < 0.65 the distribu¬ 
tion of fc -R e corrected absolute magnitudes is fairly flat. 
At redshifts z < 0.15 we have a tail of fainter objects in 
our catalogues, which are likely different types of galaxies. 
Therefore, we exclude them from this analysis. At higher 
redshifts, we start loosing fainter objects due to incomplete¬ 
ness. Since we determine both the mean luminosity and the 
mean halo mass for a given lens sample, this should not bias 
the overall mass-to-luminosity relation. 


2.2. Leasing measurement 

The shapes of the background galaxies are measured on 
images from the RGS2. Details on the data reduction 
and the shape measu rement process can be found in 
Ivan Uitert et al.l (|201lf) , and some importan t improvements 
of our lensing analysis were discussed in ICacciato et al.l 
(|2014[1 . It suffices to say that we me asure the shapes of 
the galaxies with the KSB method (iKaiser et al.l 119951 : 
iLuppino fc KaiseJ 119971 : iHoekstra et al.l 119^1. usi ng th e 
implementation described bv IHoekstra et 311 1 199811200011 . 
This method was tested on a range of simulations as part 
of the Shear Testing Programme (STEP) 1 and 2 (the 
‘HH’ method in Heymans et al. 2006 and Massey et al. 
2007 respectively) where it was found to have a multi¬ 
plicative bias of a few per cent at most and a negligi- 
ble additive bias. Recently, IHoekstra et al.l (1201511 found 
that these results were driven by the overly simplistic na¬ 
ture of the STEP simulations; for more realistic simula¬ 
tions, KSB suffers from noise bias (iKacprzak et HI 120121 : 
iMelchior fc Violall2012l : iR.efregier et al.ll2012ll as any other 
shape measurement method that is currently in use. We 
calibrate our KSB implementation on realistic image sim¬ 
ulations generated with GalSiirQ (|R,owe et al.ll2014ll that 
are set up to closely match RGS2 observations, i.e. with 
an intrinsic ellipticity distribution that is matched to the 
observations, a realistic range of Sersic prohle indices for 
the simulated galaxies, and up to a magnitude limit that 
is matched to the RCS2. We determine the multiplicative 
bias as a function of seeing: 


rucorr = -0.065 X (FWHM - 0.7) - 0.123, 


( 1 ) 


with FWHM the size of stars in an image. We use this to 
correct the shear measured in each RGS2 image. We do 
not need to correct for residual additive bias as that gen¬ 
erally averages out in galaxy-galaxy lensing due to sym¬ 
metry in lens-source pair orientations. The (multiplicative) 

^ https://github.com/GalSim-developers/GalSim 
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Redshift 

Fig. 3. The distribution of absolute magnitudes and red- 
shifts of our LRG lenses after the fc-correction (top) and 
after the /c-correction and luminosity evolution correction 
(bottom). The solid cyan line indicates the median. The 
green dashed boxes show the redshift and luminosity cuts 
on our lens sample; the magenta pentagons indicate the 
mean redshift and luminosity of those bins. The total den¬ 
sity of LRGs as a function of redshift is shown by the black 
and red histograms on the horizontal axis for the LOWZ 
and GMASS galaxies, respectively. 


noise bias correction increases the average leasing signal by 
10-15%. Note that the correction is not very sensitive to 
the adopted width of the intrinsic ellipticity distribution, 
but that it is critically important to simulate galaxies up 
to ^1.5 magnitudes d eeper than the nomin al magnitude 
limit of the survey isee lHoekstra et aT]l2015ll . 

The lensing signal is extracted by azimuthally averaging 
the tangential projections of the ellipticities of the source 
galaxies in concentric radial bins, i.e. by measuring the tan¬ 


gential shear as a function of projected separation: 




AE{R) 

L^crit 


(2) 


where AS(i?) = fi(< R) — f^(R) is the difference between 
the mean projected surface density inside radius R and the 
projected surface density at i?, and Scrit is the critical sur¬ 
face density: 

with Di, Ds, and D\s the angular diameter distance to 
the lens, the source, and between the lens and the source, 
respectively. All galaxies with an apparent magnitude of 
22 < r' < 24 and a well-defined shape measurement are 
selected as sources. 

We measure the lensing signal from the BOSS lenses in 
each 1x1 deg^ RCS2 pointing, including the sources from 
the neighbouring pointings (if present). We bootstrap over 
these patches to obtain the covariance matrix, which ac¬ 
counts for intrinsic shape noise, measurement noise, as well 
as the contribution from large-scale structures. The off- 
diagonal elements are consistent with zero on the radial 
range of interest (< 2/i^g^Mpc), hence we only use the in¬ 
verse of the diagonal as the errors on the measurement when 
fitting the models to the data. 

As these patches overlap, the contribution from large- 
scale structures might be somewhat underestimated at 
large scales; therefore, as a test, we also perform the lensing 
measurements on 2x2 deg^ non-overlapping patches and 
use that in the bootstrap resampling. We find that for 
both methods, the signal and the covariance matrix barely 
change in the radial range that we use in this work. Since 
the total area decreases if we limit our analysis to 2 x 2 deg^ 
non-overlapping patches only, and since it makes no dif¬ 
ference to the signal and its error, we decide to use all 
1x1 deg^ RGS2 pointings plus neighbours as basis for the 
bootstrapping. 

In order to compute Scrit we need the distances to the 
lenses and sources. We compute D\ for each lens separately 
using its spectroscopic redshift from SDSS. The lensing ef¬ 
ficiencies {D\s/Ds) are determined by averaging over the 
source redshift distribution, which is obtained by apply¬ 
ing the same r'-band selection to the photometric r edshif t 
catalogues of the GOSMOS field from lllbert et al.l l|201.‘lll . 
T he procedure is descri bed in more detail in Appendix C 
of ICacciato et al.l (l2014ll . Note that we previously used the 
photometric redshift catalogues from lllbert et al.l (|2009ll as 
the former was not yet publicly available. Comparing the 
average lensing efficiencies from the two catalogues, we find 
them to agree for low lens redshifts, but to be increasingly 
different at higher redshifts (up to 15% at zi=0.7). At in¬ 
creasingly high lens redshift the lensing efficiencies are more 
sensitive to the form of the adopted source redshift distri¬ 
bution, which is somewhat different for the two catalogues. 
We discuss the robustness of the derived lensing efficiencies 
in more detail in Appendix B. 

A fraction of the sources is physically associated to the 
lens galaxies, representing an overdensity of sources galax¬ 
ies that are not lensed. We cannot remove them since we 
lack redshifts for our sources. Such a contamination in the 
source catalogue dilutes the lensing signal. This can be cor¬ 
rected for by measuring the excess source number density 
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relative to the background as a function of projected sepa- 
ration, and boosting the leasing signal with this fa ctor (e.g. 
iMandelbaum et aLll2006bl : Ivan Uitert et alJ [201111 . We fol¬ 
low the same procedure. 

This boost correction itself is biased low as the galax¬ 
ies associated to the lens, and the lens itself, block light 
from the background sky, suppressi ng the source number 
densit y. The effect is described in ISimet fc MandelbaumI 
(|2014ll . As discussed in that work, a correction for this bias 
is obtained by multiplying the boost correction with a fac¬ 
tor 1/(1 — /obsc), where /obsc is the fraction of the sky that 
is obscured by the foreground galaxies. We compute this 
by using the ISOAREAJMAGE keyword in SExtractor, 
which stores how many pixels a galaxy spans on the sky. 
/obsc is taken to be the sum of these values of all galaxies 
whose centroids fall inside a radial bin, divided by the total 
number of pixels in that bin (accounting for the effect of the 
survey masks and geometry). Before doing the correction, 
we subtract the average sky-background of /obsc from the 
one observed around the lenses, as we are only interested 
in the additional obscuration. We find that the correction 
is at most 5%, in the radial bins closest to the most lumi¬ 
nous and low-redshift lenses. The correction decreases at 
larger radii, as well as for fainter, higher redshift LRGs, as 
expected. 

The robustness of the leasing s i gnal h as been addressed 
in Appendix B of ICacciato et all (l2014ll . There, we show 
that the cross shear signal of our lens sample is consistent 
with zero. Also, the random shear signal, which is used to 
correct for the effect of residual systematics in the shape 
measurement catalogues, is smaller than the real signal for 
all the projected separations we use in this work. However, 
an overall multiplicative bias could still be present, either 
through an incorrect determination of the noise bias cor¬ 
rection, or through the use of incorrect leasing efficiencies. 
In Appendix B we perform an internal consistency check of 
our measurement pipeline, which provides strong evidence 
that such a bias is unlikely to be significant. 


3. Halo model 


In this section we describe the model that we employ 
to provide a physical interpretation of our measurements. 
The halo model provides a useful framework to describe 
the stacked weak leasing signal around galaxie s (see 
e.g. Maiidglbaum_e1^aLl l2006bl : iGacciato et al. l2009l I20l4 


iMivatake et al.ll20I,lll . It is based on a statistical descrip¬ 
tion of dark matter properties, such as their average density 
profile, their abundance, and their large scale bias, comple¬ 
mented with a statistical description of the way galaxies of 
a given luminosity populate dark matter haloes of differ¬ 
ent masses (also known as halo occupation statistics). In 
its fundamental assumptions , the model is similar to the 
one presented i i ilSeliakI (|200Clfl , ICoorav &: ShethI (|2002ll and 
IGacciato et al.l (|2009fl . 

Galaxy-galaxy leasing probes the average matter distri¬ 
bution projected along the line-of-sight at a given projected 
physical separation, i?, for a set of lenses. The quantity of 
interest is the excess surface mass density profile, AE(i?), 
which is determined from the projected surface mass den¬ 
sity, T,{R). Since we measure the average signal of many 
lenses, the projected matter density can be expressed in 


terms of the galaxy-dark matter cross-correlation, ^gm('c): 

m=Pr. [l+4„.(r)]db;, (4) 

Jo 

where the integral is along the line of sight, w is the comov¬ 
ing distance from the observer, Wg the comoving distance to 
the source and is the mean matter density at the red- 
shift of the lens. The three-dimensional comoving distance 
r is related to w through = ta/ + uj'^ — 2 ujiujcos9, with 
wi the comoving distance to the lens and 9 = R/D\ the an- 
gular separation betw een lens and source (see Figure I in 
ICacciato et al.l (1200911 b Note that the galaxy-dark matter 
cross-correlation is evaluated at the average redshift of the 
lens galaxies. 

Under the assumption that each galaxy resides in a dark 
matter halo, AE can be computed using a statistical de¬ 
scription of how galaxies are distrib uted over dark mat¬ 
ter h aloes of different masses (see e.g. Ivan den Bosch et'HI 
1201,41 . Specifically, it is fairly straightforward to obtain the 
two-point correlation function, Cgm(?', ,z), by Fourier trans¬ 
forming the galaxy-dark matter power-spectrum, Pg^(k, z), 
i.e. 

= Pgm{k,z)^^^^k'^dk, (5) 

with k the wavenumber. The quantity Pgm(fc, z) can be ex¬ 
pressed as a sum of a term that describes the small scales 
(one-halo, Ih) and one that describes the large scales (two- 
halo, 2h), each of which can be further subdivided based 
upon the type of the galaxies (central or satellite) that con¬ 
tribute to the power spectrum. This reads 

P^mik) = P//)(fc) + P//((fc) + P2h(fc) + p2^(fc) . (6) 

The terms in Equation [5] can be written in compact form 
as 

p/y(/c,2;)= J 'Hx{k,M,z)'Hyik,M,z)nh{M,z)dM, (7) 

P^^{k,z) = J dMiR^^k, Ml, z) ni,{Mi,z) 

J dM2Py(fc,M2,z)nh(M2,z)Q(fc|Mi,M2,z), (8) 

where ‘x’ and ‘y’ are either ‘c’ (for cen¬ 
tral), ‘s’ (for satellite), or ‘m’ (for matter), 
Q{k\Mi,M 2 ,z) =bh{Mi,z)bh{M 2 ,z)P^(k,z) describes 
the power spectrum of haloes of mass Mi and M 2 , 
and i t contains t h e larg e scale bi as of haloes bh(M) 
from IrToker et all (|2010ll (but see Ivan den Bosch et"^ 
(I 2 OI 3 D for a more sophisticated mo delling of this term) . 
ni,{M, z) is the halo mass function of iTinker et al.l i|2010ll . 
Furthermore, we have defined 


M 

Rmik, M, z) = — Uh(fc|M, z), 

(9) 

Pm 


nAk,M,z)= Ux(fc|M) 

n^{z) 

(10) 


where 

Uc{k\M) = 1-poff+Poff xexp[-0.5fc^(rs{M}7^off)^], (H) 
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and 

Us{k\M,z) = u\i{k\M,z) . (12) 

PoS is the parameter that describes the probability that 
the ‘central’ galaxy does not reside at the centre of the 
dark matter halo, whereas T^off quantifies the amount of 
off-c entering in terms of the halo scale ra dius, rs{M) (see 
e.g. ISkibba et alJ 120111 : iMore et ^ l2014ll . In our fiducial 
model, we set Poff = = 0 , but we explore the impact 

of this assumption in Section 0] The functions {Nc\M) and 
{Ns\M) represent the average number of central and satel¬ 
lite galaxies in a halo of mass M = 47r(200pm)^20o/3j 
fined as: 


1 / (logM-logMmean)^ 

{Nc\M) = ,— -exp-—o- 

V^ln(10)McriogM \ 

(13) 

(iVs|M) = (M/Ml) /tran.(M) (14) 

where 


/trans(M) = 0.5 X 


1 -I- erf 


logM - log Meat) 

^ trans 


• (15) 


We use a flat, non-informative prior for uiogM and Mmean, 
set Meut = (Mefi) (see Equation [TH]) and Utrans = 
0.25. Since LRGs ar e thought to be predominantly cen¬ 
tral galaxies (see e .g. IWake et "ni 120081 : IZheng et al.l 120091: 
iPareiko et al.ll20i^ . we set (A^s|M) to zero and only fit for 
Mmean and (Tm in our nominal runs. We test the impact 
of this assumption on the derived quantities in the result 
sections by additionally fitting for Mi. We have tested that 
the result is fairly insensitive to the details of the modelling 

of /trans(M). 

hg(z) is the number density of galaxies at redshift z: 


ng(z) = J (A^g|M)nh(M, z)dM 

« /(iVc|M)nh(M,z)dM. 


(16) 


The last equality is exact in the case of LRGs being only 
central galaxies. uii{k\M, z) is the Fourier transform of the 
normalized density distribution of dark matter within a 
halo of mass M, for which we assume a N avarro-Frenk- 
White (NFW) profile (iNavarro et al. 199611 a nd a mass- 
concentration relation from Duffy et al.l (|200^ : 


c =a(^ 
“ VA^pivot 


(1 + ^)'', 


(17) 


with A = /cone X 10.14, B = —0.081, C = —1.01, and 
Mpivot = 2 X lO^^h“^M 0 . Note that /cone is a free param¬ 
eter that allows the normalisation of this relation to vary. 
Specifically, we apply a non-informative flat prior on this 
parameter. 

The average halo mass in a given luminosity bin, which 
is what we shall refer to as ‘effective’ halo mass in what fol¬ 
lows, can then be computing taking into account the weight 
of the halo mass function: 


, /(A^c|M')n^(M',^iens)M'dM' 

^ /(iVc|M')n,,(M',Ziens)clM' 


(18) 


where ziens is the mean redshift of the lens galaxies in a 
given luminosity bin. The distinction between the mass as¬ 
sociated with the mean of the log-normal distribution and 


the halo mass inferred accounting for the mass function is 
of relevance because LRGs populate fairly massive haloes 
for which the mass fun ction is steep (see e.g. Figure 7 in 
iLeanthand et al.l (1201511 1. 

At small scales one expects the baryonic mass of LRGs 
to contribute to the lensing signal. The smallest scale used 
in this study is 50 h/g^kpe, which is much larger than 
the typical extent of the baryonic content of a galaxy. 
Therefore, it is adequate to model the lensing signal of the 
LRGs itself as a point source of mass Mg « M*. This reads 

AEih.g(i?) ^ . (19) 

) We use the value of the average stellar masses, (M*), for 
the galaxies in the luminosity bins under investigation here. 
The stellar masses are obtained by matching our lens cat¬ 
alogue to the MPA-JHU stellar mass catalogue^- As the 
MPA-JHU catalogue is based on the MAIN sample from 
SDSS, we only have matches at low redshift. However, the 
point mass has a small impact on our fit results; hence we do 
not expect that a potential evolution of the average stellar 
mass-to-light ratio for LRGs can be so strong that it could 
significantly affect our results. 

To summarize, our fiducial model for the lensing sig¬ 
nal is the sum of three terms: one describing the lensing 
due to the baryonic mass; the second is responsible for the 
small (sub-Mpc) scale signal mostly due to the dark mat¬ 
ter density profile of haloes hosting central LRGs; and the 
last describes the large (a few Mpe) scale signal due to the 
clustering of dark matter haloes around LRGs. This reads: 

AE(i?) = AEi''’S(R) + + AS^^ (A). (20) 

We simultaneously fit the halo model to the three luminos¬ 
ity bins and do this for each redshift slice separately. We 
have five free parameters in each fit: the three mean masses 
of the luminosity bins, the scatter and the normalisation of 
the mass-concentration relation. The fit is performed using 
a Markov Chain Monte Carlo Method (MCMC). Details of 
its implementation can be found in Appendix [Cl 

We fit the model to the measurements on scales be¬ 
tween 0.05 and 2 Mpc. At these scales, both the mea¬ 
sured lensing signal and the halo model predictions are 
fairly robust. At larger scales, the lensing signal becomes 
smaller and more susceptible to residual systematics. At 
scales smaller than 0.05 h/g^Mpe, lens light may bias the 
shape measurements. For the halo model, the overlap be¬ 
tween the 1 -halo and 2 -halo term is notoriously hard to 
model because of, amongst others, halo exclusion and non¬ 
linear biasing. This mainly impacts the few Mpc regime. 
Most of the information about the halo masses and con¬ 
centrations is contained in the lensing signal within the 
virial radius, so we do not loose much statistical precision 
by limiting ourselves to these scales. 

4. Luminosity-to-halo mass relation 

To study how the luminosity-to-halo mass relation of LRGs 
evolves with redshift, we divide our sample in bins of (k+e)- 
corrected absolute magnitude and redshift as detailed in 
Tableland Figure [S] For z < 0.43, we only select lenses 
from the LOWZ sample, at higher redshifts we exclusively 

® http://www.mpa-garching.mpg.de/SDSS/DR7/ 
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Fig. 4. The lensing signal AS of LOWZ (top two rows) and CMASS (bottom two rows) lenses as a function of projected 
separation for the three luminosity bins (after the (fc + e)-correction is applied). The solid red lines show the best-fit halo 
model, the orange and yellow regions the 1 and 2a model uncertainty, respectively. We fit the signal on scales between 
0.05 and 2 Mpc. 


select CMASS galaxies. The average log stellar masses for 
the consecutive luminosity bins are 11.2, 11.5 and 11.7 
[\og{h^Q Mq)]. For each bin we measure the average lensing 
signal, which is shown in Figure |4j together with the best- 
fit halo models and the model uncertainties (computed as 
detailed in Appendix ICl) . We find Xred values of 1.8, 1.6, 
1.2 and 1.0, going from the lowest to the highest redshift 
slice. Hence the fits of the CMASS samples are good, but 
for the LOWZ samples the Xred values are somewhat large, 
suggesting that either our error bars are underestimated, or 
that the model that we ht to the data is overly simplistic. 

The errors on the lensing measurements account for 
intrinsic shape noise, measurement noise and the impact 
of large-scale structures. We have, however, ignored some 
small sources of error, as their amplitude is much smaller 
than the statistical errors on the lensing signal: the error 
on the boost correction, which is typically a few percent at 
small scales; the error on the obscuration correction, which 
is even smaller; the error on determining the lensing effi¬ 
ciency and the error on the multiplicative bias calibration, 
whose magnitudes are unknown but are likely of the order 
of a few percent. Combined, they might increase the errors 
by as much as ~10%, although the exact number is diffi¬ 
cult to estimate reliably. If we increase our error bars by 
this amount, we would get Xred values of 1.5, 1.4, 1.0 and 
0.8, respectively. The fact, however, that we find reasonable 
Xred values for the CMASS sample, but not for the LOWZ 
sample, suggests that a systematic underestimate of our er¬ 
rors is unlikely to be the dominant cause. 

Even though a visual inspection of the covariance 
matrix lead us to believe that it is diagonal on scales 
<2h^Q Mpc, there could be low-level off-diagonal terms 
present that, if included, would lower the Xred values. This 


potentially affects the LOWZ results more, as the measure¬ 
ments have a higher signal-to-noise ratio and the covariance 
matrix is less noisy. To test this, we recompute the Xred val¬ 
ues using the full covariance matrix which we obtained from 
bootstrapping (see Section l^^ for the best-fit models. Note 
that we only include the covariance between radial bins of a 
lens sample, but not the covariance between the radial bins 
of the different luminosity samples. If present, they would 
lower the Xred values even more. We Hnd that Xred 
first redshift slice reduces to 1.6, while it does not change 
for the other three redshift slices. Hence the effect is small 
and does not fully explain the high Xred values. 

Figure 0] shows that the signal-to-noise ratio of the lens¬ 
ing measurements of the LOWZ samples is very high and 
would allow for a more sophisticated modelling. When we 
include a satellite term or a miscentering term in the halo 
model, however, the Xred values do not improve, as the 
lensing signal alone cannot constrain the miscentering dis¬ 
tribution parameters very well, and the expected number 
of satellites is low. Allowing for even more freedom in the 
fit might lead to overfitting of the CMASS results. Using 
different halo models for the different samples, or splitting 
the LOWZ sample up in more luminosity bins, reduces the 
homogeneity of the analysis, which is one of the key advan¬ 
tages of our work. Hence we choose to stick to the settings 
described above. In Section lO we perform a sensitivity 
analysis and find that our results do not critically depend 
on various choices in the analysis, suggesting that the quan¬ 
tities we derive from the fits are robust. 

In the halo model, we fit the mean and the scatter of 
the log-normal distribution that describes {Nc\M). The log 
of the mean has typical values of ~14.5, ^15 and ^15.5 for 
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Table 1. Properties of the lens bins (after {k + e)- 
correction) 



Nlens 

(z) 

{Lr) 

Meff 

fconc 

Xred 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

0.15<2<0.29 (LOWZ) 






Llzl [-21.8,-21.2] 

2969 

0.219 

0.65 

3-65«jg 



L2zl [-22.4,-21.8] 

2606 

0.226 

0.99 

5.60l°;S 

0.681^6 

1.8 

L3zl [-22.8,-22.4] 

300 

0.234 

1.58 

14.9tl:l 



0.29<z<0.43 (LOWZ) 






Llz2 [-21.8,-21.2] 

3771 

0.351 

0.66 

^•oO_Q_45 



L2z2 [-22.4,-21.8] 

4502 

0.364 

1.00 

c 1 r:r+0.69 
>J--LO_o.64 

0.85l°:l? 

1.6 

L3z2 [-22.8,-22.4] 

721 

0.368 

1.58 

9.86ll;«® 



0.43<2:<0.55 (CMASS) 



2.03t°:« 



Llz3 [-21.8,-21.2] 

8530 

0.499 

0.61 



L2z3 [-22.4,-21.8] 

4213 

0.503 

0.99 

_Q YY 

0.77t°:“ 

1.2 

L3z3 [-22.8,-22.4] 

587 

0.500 

1.59 




0.55<2<0.70 (CMASS) 



1 Q9+0.66 

i.yz_Q 57 



Llz4 [-21.8,-21.2] 

5256 

0.596 

0.64 



L2z4 [-22.4,-21.8] 

5161 

0.611 

1.00 

4.16tJ;°J 

n 'To+u.i^io 
' '^-0.20 

1.0 

L3z4 [-22.8,-22.4] 

969 

0.616 

1.60 

6.48t?:“ 




Notes. (1) absolute magnitude range (after (fc + e)-correction); 
(2) number of lenses; (3) mean redshift; (4) mean luminos¬ 
ity [IO^^/iyq^Lq] (after (fc -|- e)-correction); (5) best-fit halo 
mass [10^^Mq]; (6) best-fit normalisation of the mass- 
concentration relation; (7) reduced chi-squared of the fit. 


the three luminosity bins, while the scatter ranges between 
0.7 and 0.8. Neither evolves with redshift. Note that this is 
the scatter in the log of the halo mass and no t in luminos¬ 
ity. The latter was fit in ICacciato et al.l (1201411 where it was 
found to have a value of CTiogLo = 0.146 ± 0.011, obtained by 
fitting the halo model to the lensing signal of all galaxies 
in the DR9 that overlap with RCS2. The scatter in halo 
mass is much larger, because the luminosity-halo mass 
relation flattens at higher luminosities; a small scatter in 
luminosity corresponds to a large one in halo mas s (see 
e.g. Figure 3 and the discussion in iMore et al.l (l2009ll l. 

The quantity of interest that we can compare to other 
works is the ‘effective’ halo mass, which is given in Table 
[TJ We plot it as a function of luminosity in Figure [S] 
We Hnd that the masses increase with luminosity and 
decrease with redshift. To quantify this, we parametrize the 
luminosity-to-halo mass relation by M^e = Mo^l{L/Lo)^’^ , 
using a pivot luminosity of Lq = The best-ht 

power law fits are shown in the same figure; the confidence 
contours of the fitted amplitude and slope are shown in 
Figure El We have listed the fit parameters in Tabled 
The slope of the luininosity-to-halo mass relation does 
not change significantly for our different redshift samples 
and has a typical value of 1.4. The amplitude, however, is 
about ~4tT higher for our lowest redshift slice compared 
to the highest one. On average, the masses of LOWZ 
galaxies increase by 251^4% between redshift 0.36 and 
0.22; the masses of CMASS galaxies increase by lOl^o^o 
from redshift 0.6 to 0.5. If we assume that CMASS galaxies 



Fig. 5. The mean (fc -|- e)-corrected luminosity versus the 
best-fit halo mass of the twelve lens samples. The different 
symbols indicate the different redshifts bins, as indicated in 
the figure. The coloured areas indicate the 68% confidence 
intervals of the powerlaw fits. 



^eff(%=10" L^) 

Fig. 6. 68% confidence contours for two parameters of the 
powerlaw fit between luminosity and halo mass. 


evolve into LOWZ galaxies and combine the results, we 
hnd an average increase of 80 ^ 28 ^ (stat. errors) in Mes 
at Lq = IO^^H^qLq from z^O.6 to z~0.2. Fixing the slope 
to its average value of 1.4 and only htting the amplitude 
ch anges this number to 77^27^- 

iToieiro et al.l (|2012f) found that at brighter absolute 
magnitudes, a larger fraction of CMASS galaxies are the 
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progenitors of LOWZ galaxies. If we discard the lowest 
luminosity bin, we find that, for a pivot luminosity of 
Lq = 1-3 X Lq, the average halo mass increases 

with 97^38%- Considering only the brightest luminosity 
bin, the average halo mass increases with 160 ^ 7 g^%. 


4.1. Sensitivity analysis 


To study how sensitive our results are to the adopted 
luminosity evolution correction, we perform the following 
test. We multiply our nominal correction with a factor, 
such that the resulting luminosity evolution curves roughly 
cover the range of reasonable models that are shown in 
the lower panel of Figure O The factors we choose are 
1 — 0.2 X 2 and 1 + 0.2 x z, respectively. Next, we recompute 
the luminosities, repeat the lensing measurements (using 
the same cuts) and the halo model fits, and compare the 
resulting best-ht halo masses. We find that the best-ht 
halo masses of the individual luminosity bins do not signif¬ 
icantly shift compared to the nominal results. We fit the 
luminosity-to-halo mass relation and list the parameters in 
Table [2] As is shown there, we find that the best-fit slopes 
are very similar. The best-fit amplitude shifts with Icr at 
most compared to our nominal results. For the 1 — 0.2 x z 
modification factor, the resulting increase in halo mass of 
LRGs from z = 0.6 to z = 0.2 is 85^30%; for the 1 -|- 0.2 x z 
modification factor, it is SOl^s^o- The halo masses and 
corresponding growth change somewhat because the lens 
selection shifts systematically, such that we analyse lens 
samples that to some extent are intrinsically different. 
Importantly, however, our results do not critically depend 
on the choice of the luminosity evolution correction. For 
future work that is expected to have an improved statistical 
precision, more detailed knowledge of this correction will 
be required. 

To test how sensitive our results are to the assumption 
that all LRGs are located at the centre of their dark matter 
haloes, we do two halo model runs where we allow for more 
flexibility. First, we assume that a fraction of the LRGs is 
miscentered, following Equation[TT] We use PoS and i?off as 
additional free parameters with a flat uninformative prior 
in the range [0,1]. The allowed miscentering distribution 
ranges from the lenses being all correctly centered (poff 
=i?off = 0 ) to all being miscentered and located at the 
halo scale radius (poff =Roe =!)■ We find PoS = O.STIqJ®, 
0 . 26 ^ 0 , 23 ) 0-50lQ 3 g and O.dOtg^go for the low and high 
redshift bins of LOWZ and GMASS, respectively, with 
corresponding miscentering radii of i?off = 0 . 371 q 25 , 

0.32^0 29 , 0 . 491 qJ 7 and The resulting powerlaw 

parameters are listed in Table [2] We find that the powerlaw 
parameters do not change significantly. The total increase 
in halo mass of LRGs corresponds to 75227 %: consistent 
with our nominal result of 80228%- These constraints on 
the miscentering distribution are in broad agreement with 
previous galaxy-galaxy lensing and clustering result s of 
GMASS galaxies. For instance, iMivatake et al.l ( 2013f) re¬ 
port PoS < 0.66 and i?off = 0 . 792 o 38 : whereas [More et ^ 
(12014 find PoS = 0.34 ± 0.18 and RoS = 2.22];^. 

As a related test, we study how our results change 
when we assume that a fraction of LRGs are satellites. We 
use a simple model with only one free parameter, i.e. Mi. 


Table 2. Powerlaw parameters, normalisation of the mass- 
concentration relations and reduced chi-squared values for 
various runs as described in the text. 


Mo.l 

nominal run 

0.15< z <0.29 6.43 ±0.52 

0.29< z <0.43 5.14 ±0.44 

0.43< z <0.55 4.08 ±0.51 

0.55< z <0.70 3.74 ±0.60 

e-correction x (1 — 0.2z) 
0.15< z <0.29 5.90 ±0.50 

0.29< z <0.43 5.03 ± 0.45 

0.43< z <0.55 3.82 ± 0.47 

0.55< z <0.70 3.33 ±0.58 

e-correction x (1 ± 0.2z) 
0.15< z <0.29 6.12 ±0.51 

0.29< z <0.43 5.51 ±0.47 

0.43< z <0.55 4.59 ± 0.56 

0.55< z <0.70 4.15 ±0.68 

miscentering run 
0.15< z <0.29 6.10 ±0.50 

0.29< z <0.43 5.08 ± 0.44 

0.43< z <0.55 3.87 ±0.49 

0.55< z <0.70 3.75 ±0.60 

satellite fraction run 
0.15< z <0.29 4.93 ±0.45 

0.29< z <0.43 4.22 ± 0.38 

0.43< z <0.55 3.17 ±0.35 

0.55< z <0.70 2.88 ±0.44 


Pl 

fconc 

Xred 

1.59 ±0.29 

0 . 6822 o 6 

1.8 

1.42 ±0.27 

0.8522.n 

1.6 

1.27 ±0.31 

0.772215 

1.2 

1.32 ±0.44 

'70+0-25 
' '^-0.20 

1.0 


1.56 ±0.29 

0.7422J? 

1.9 

1.44 ±0.27 

0.8322m 

1.3 

1.19 ±0.35 

0.8022m 

1.0 

1.43 ±0.53 

0.75222° 

1.0 

1.65 ±0.30 

0.7522J? 

1.9 

1.53 ±0.26 

0.8422m 

1.4 

1.35 ±0.31 

0.8322m 

1.3 

1.58 ±0.49 

r\ '70 + 0.25 
' '^-0.20 

1.0 

1.61 ±0.29 

0.922211 

1.9 

1.41 ±0.27 

0-9522m 

1.8 

1.26 ±0.31 

1.04221? 

1.3 

1.31 ±0.44 

0.932218 

1.1 

1.39 ±0.26 

0.8722m 

1.9 

1.27 ±0.25 

-1 1 q+0.18 

-*-•-‘-^-0.17 

1.7 

1.02 ±0.29 

0.902211 

1.2 

1.08 ±0.43 

0.822211 

1.0 


We follow this approach because we are not interested to 
determine the satellite HOD (as lensing alone is not very 
sensitive to that), but because we want to get an estimate 
how much our results could be affected by ignoring the 
contribution of satellites. The typical satellite fraction 
for LRGs is ~10% and decreases for more massive LRG s 
(I White et al.l 1201 ll : iPareiko et al.l 1201^ iMore et ^I20l4 . 
We therefore put a prior on Mi such that the resulting 
satellite fractions from our model are between 5% and 15%. 
The resulting luminosity-to-halo mass relation parameters 
are listed in Table |2l The best-fit slopes are consistent 
with our nominal run, but the amplitude decreases by 
1-2 ct. The total increase of halo mass is 822^9% over the 
full redshift range of our LRG sample, consistent with 
our nominal result. The normalisation decreases because 
the satellites are associated to more massive haloes, with 
correspondingly larger lensing signals. This lowers the 
required contribution to the total signal from central 
LRGs, and hence their mass. This does not happen in 
the miscentering run, where the lensing signal is merely 
smeared out, but the integrated signal (and hence the 
mass) stays the same. 
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4.2. Comparison to previous work 

Several previous works have provided mass estimates of 
LRGs, using gravitational lensing, clustering, abundance 
matching, or a combination of these. The selection of the 
samples, the models fit to the data and the definitions 
of mass, generally differ between these studies, limiting 
the level of detail with which we can perform a comparison. 


4.2.1. Lensing results 

iMandelbaum et all (|2006all measured the masses for a 
sample of over 4 x 10"^ LRGs with spectroscopic redshifts 
from SDSS-I/II using weak lensing. Bluer and fainter 
LRGs are discarded using colour-magnitude cuts, as 
well as LRGs that are likely satellites of larger systems, 
hence the selection is not identical to ours. The resulting 
sample is split into a faint and bright part using a cut 
at Mr = —22.3, and the mean luminosity of the two 
samples is 5.2 x 10^° and 8.6 x lO^°h“^L 0 , respectively. 
The luminosities are computed in a similar manner as 
in our work. All LRGs are selected in the redshift range 
0.15 < z < 0.35 and have a mean effective redshift of 
0.24. Masses are estimated using NFW fits plus a baryonic 
component, where the fitting range is restricted to small 
scales were the 2-halo term can be neglected. The masses 
are defined as the enclosed mass in a sphere where the 
density is 180 times the mean background density, Migofe, 
instead of 200 times the mean density which we use; the 
difference between these definitions is only a few percent. 
The faint LRGs are found to reside in haloes of masses 
2.9 ± 0.4 X lO^^/i“^M0 and the bright ones in haloes 
of masses 6.7 ± 0.8 x lO^^/i“^M 0 . The measurements 
are shown in Figure El We find that the masses of our 
low-r edshift slices are su b stantia lly larger than the ones 
from IMandelbaum et ^ (|2006ari . The discrepancy may 
be caused by the scatter be tween luminosity and hal o 
mass, which is not included in IMandelbaum et al.l (|2006af) . 
Allowing for a non-zero scatter results in larger halo 
masseS; _ 

In iMivatake et al.l (|2013ll . the lensing signal of 4,807 
CMASS galaxies with 0.47 < z < 0.59 is measured in 
the overlap with CFH TLS using the publ icly available 
CFHTLenS catalogues (|Hevmans et al.l[2012ll . The lensing 
signal is ht together with the projected clustering signal us¬ 
ing a halo model that is similar to the one we have adopted 
here. Halo masses are defined with respect to 200 times 
the background density as we do. The average halo mass 
of GMASS galaxies is found to be 2.3 ± 0.1 x Mq. 

To compare it with our results, we compute the average 
luminosity of CMASS galaxies with 0.47 < z < 0.59 in 
our catalogue. We find a value of 3.8 x lO^°/i“^L 0 and 
assume that this number is representative for the average 
luminosity of the lenses in that work. This estimate is 
actually a bit too low, as IMivatake et al.l (|2013ll also apply 
a cut on stellar mass to remove the least massive (hence 
faintest) objects, which we cannot mimic. However, as an 
indication, we find that if we remove the faintest 10 % of 
our CMASS sample, the average luminosity only increases 
to 4.1 X lO^°/i“^L 0 , hence it is unlikely that the average 
luminosity is very far off. We compare the results in Figure 
El and find that our CMASS masses are somewhat lower 
but consistent. 


Several papers have studied how the average mass 
of galaxies changes as a function of redshift (e.g. 
van Uitert et al.l 120111 : IChoi et ^ 120121 : Leauthaud et al.l 


I20ia 

20121: iTinker et al.l 120131: iHndsou et a.l.ll2fl15ll . using weak 

lensing measurements. Since these works do not target 
LRGs specihcally, and since the modelling of the signal 
differs from our approach, we cannot compare th e re- 
sults in deta il . How ever, both iTinker et al.l l|2013l1 and 
iHudson et al.l (|201,5^ include red, massive galaxies in their 
wo rk, hence we c an at least compare the recovered trends. 

ITinker et al.l (l2013l1 use measurements of weak lensing, 
clustering and the stellar mass function of galaxies in 
COSMOS to constrain the stellar-to-halo mass relation. 
Halo masses are defined like ours. The relations they report 
predict the average logio(Af*) as a function of halo mass, 
instead of the average halo mass at a given stellar mass, 
which is what we measure; these relations are not the same 
due to intrinsic stellar mass scatter, which is illustrated 
in their Figure 7. From the right-hand panel of that 
figure we observe that for passive galaxies, at the average 
stellar masses of our samples, the mean halo mass is 
roughly ~0.5 dex lower than what we hnd. There are many 
differences between the analyses that could contribute to 
this difference, such as systematic offsets between stellar 
mass estimates, the selection of the samples, the definition 
of mass and the modelling of the signal. Nonetheless, at a 
stellar mass of logiQ(M*) ^ 11.4 (typical for LRGs), the 
average halo mass increases from ^ 10 ^^ ® to Mq 

between redshifts of 0.88 and 0.36, an increase of almost 
100 %, similar to the average growth in halo mass that we 
fin d for our LRG s ample from redshift 0.62 to 0.21. 

iHudson et al.l (|2015l 1 use the shape measurements 
from CFHTLenS to measure the lensing signal for blue 
and red galaxies in three redshift slices. The lenses are 
binned in luminosity rather than stellar mass to avoid an 
Eddington bias due to the larger observational errors on 
stellar mass compared to luminosity. The stellar mass is 
then determined using the mean stellar-mass-to-luminosity 
ratio. For the highest luminosity bin of red lenses, which 
has a mean stellar mass of ~2 x Mq, the average 

halo mass increases from 0.84 ± 0.17 x Mq to 

1.32 ± 0.33 X W^H^^Mq from redshift 0.67 to 0.29. 
The masses are defined with respect to pcrit instead of 
the mean density, resulting in masses that are ^30-40% 
smaller compared to ours. Furthermore, the intrinsic 
scatter between luminosity/stellar mass and halo mass is 
not accounted for in their modelling, which also leads to 
smaller masses. Finally, the selection of the lens samples 
differs. However, the ~58% increase in average halo mass 
is comparable to what we find. 


4.2.2. Clustering results 


The clustering of LRGs is well-studied in the literature and 
has been used to derive halo masses (e.g. Blake et al. 2008: 


Wake et al.ll2008HZheng et al.ll2009HSawanefwit et al.ll 2011 


Nikoloudakis et ahl l2013l: Pareiko et al.l l2013t iGuo et al.l 
201411 . iPareiko et al.l ( 201311 measure the clustering of 
galaxies with 0.2 < z < 0.4 from the LOWZ sample and 
fit it with a halo model. The probability distribution of 
halo masses, as shown in their Figure 9, has a mean of 
5.2 X lO^^/i“^M0. We plot it in Figure El and find that 


11 














































































Van Uitert, Cacciato, Hoekstra & Herbonnet: Luminosity-to-halo mass relation of LRGs 


it is larger than our LOWZ measurements. Although not 
specified, we assume that the mass is defined as Misnt, a s 
is mentioned in a companion paper (|White et al.l l201lli , 
which is comparable to our definition. Their halo mass 
distribution is fairly broad, however, and our constraints 
ma y well fall i n side t heir 68 % confidence region. 

IGuo et id] (|2014ll measure the clustering of CMASS 
galaxies divided in three i-band magnitude selected 
samples, which we cannot directly compare to our mea¬ 
surements. They find that going from their faintest to 
their brightest bins, the peak host halo mass increases 
from 1.1 X to 3.3 x which is 

quite comparable to our results. Their masses are defined 
with respect to the mean density, like ours. However, we 
measure an ‘effective’ mass and not the peak halo mass, 
an d it is unclear how much these definitions differ. 

IZheng et al.l (|2009ll fit the clustering signal of SDSS-1/11 
LRGs using an HOD approach. Their LRGs are divided 
in a faint and bright sample using their Mg magnitude, 
hence we cannot directly compare. The distribution of 
halo masses (defined like ours) of these samples peaks at 
~4.5 X 10^^/i“^Mq and ~lO^'^/i“^M 0 , respectively, similar 
to the values we find for faint and bright LOWZ samples, 
but note, again, that we measure the effective mass and 
not the peak halo mass. The scaling of luminosity with 
host halo mass, which is given by Lc oc is consistent 

wit h our result s _ 

I Wake et"^ (|2008ll measure the evolution of the clus¬ 
tering signal of galaxies from SPS S and the 2dF-SDS S 
LRG and QSO Survey (2SLAQ; iGannon et al.l l2006[) . 
They match the selections using colour and magnitude 
cuts, which complicates the comparison with our results. 
However, they find that the effective halo masses increase 
with ~50% from z = 0.55 t o z = 0.2, consistent with 
our findings. ISawangwit et al.l (| 2011 fl study three separate 
LRGs samples in SPSS, with a mean redshift of 0.35, 
0.55 and 0.75. After applying additional selection criteria 
such that the space density of LRGs is similar to the 
SDSS-I/H LRG sample, the effective halo mass, defined 
as the virial mass, is found to be 6.4 ± 0.5 x 10 ^®/i“^Mq, 
4.7 ± 0.2 X and 4.3 ± 0.2 x 

respectively. We assume that the average luminosity of 
these three samples is similar to that of the LOWZ sample 
and show their measurements in Figu re E] We find that 
thei r results agree fairly well w ith ours. I Blake et al.l l|2008ll 
and iNikoloudakis et all (|2013ll study the clustering signal 
of different samples of LRGs, which is difficult to compare 
with our results. The typical halo masses of LRGs are 
found to be a few times lO^®/i“^M 0 , hence broadly in 
agreement. 


4.2.3. Abundance matching results 

iMasaki et al.l (^2013^ apply an abundance matching tech¬ 
nique to N-body simulations in order to construct mock 
LRG samples whose number density matches that of LRGs 
in SPSS. From these samples a mock lensing signal is con¬ 
structed, whose shape agrees well with, but whose ampli¬ 
tude is ^ 20 % larger tha n, the measured len s ing si gnal of 
SPSS LRGs presented in iMandelbanm et al.l (l2013ll . They 
Hnd that their LRGs reside in haloes with a mean virial 
mass of 5.6 ± 0.1 x lO^®h“^M 0 . This definition of mass 
is ~ 10 % smaller than M 200 for masses and concentrations 
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Fig. 7. Evolution of the amplitude of the powerlaw fit 
between luminosity and halo mass with redshift. The 
black das hed lines shows the predicted trend from pseudo¬ 
evolution (|Diemer et al.l[20T^ for halo masses that are typ¬ 
ical for LRGs, scaled to overlap with our first data point. 


that are typical for LRGs at 2 ~0.3. We assume the aver¬ 
age luminosity equals that of SDSS-I/H LRGs and show the 
measurement in Figure [5l We find that this mass estimate 
is in fair agreement with our results. 

4.3. Interpretation 

We find that M^s of LRGs increases by approximately 80% 
from redshift 0.6 to 0.2. This is not only because of dark 
matter accretion. Part of the growth can be attri buted to 
the so-called pseudo-evolution (iPiemer et al.ll 20 l'^ . M 20 Q is 
dehned with respect to a reference density, which is redshift 
dependent. Even if a halo is static, i.e. it does not accrete 
anything, M 200 increases towards lower redshift. The halo 
mass function that we use to determine uses masses 
that are defined with respect to Pm{z). Hence the evolution 
of the halo mass function is a mix of physical evolution and 
pseudo-evolution. Therefore, the evolution of M^g is also a 
mix of the two. 

Since the halo masses of our LRGs span a relatively nar¬ 
row range, we can estimate the contribution from pseudo- 
evol ution to our observed increase in mass using the results 
from iDiemer et al.l (|2013ll . For halo masses that are typical 
for our LRGs, A /200 increases by approximately 33% from 
z = 0.6 to z = 0.2 due to pseudo-evolution, as is illustrated 
in Figure [71 If we fit the amplitude of the pseudo-evolution 
curve to our measurements, we find that Xrad = pro¬ 
viding weak support that the slope is steeper and requires 
additional dark matter accretion. 

It is interesting to compare our growth rates to those ob¬ 
tained from simulations. Ideally, we would like to compare 
to hydrodynamical simulations, but those are not available 
for the mass range we are interested in. However, a compari- 
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son with dark-matter only simulations is interesting as well, 
as a good/poor agreement points towards the relevance of 
baryonic physics. We first compare to the two Millennium 
simulatio ns, for which growth a nd merger rates have been 
derived in iFakhouri et ahl (l2010[l . They find that the growth 
rate is well described by: 

(^)mea„ = 46.1MQyr 

x(l-bl.ll2)x/fMrT^)3TT^. (21) 

Using this equation, we find that a halo of mass 
grows by 38% from redshift 0.6 to 0.22, whi lst a halo of mass 
5 X IO^^Mq grows by 43%. More recently. IWetzel &: Nagail 
(IMl use hydrodynamical simulations to measure the 
amount of physical accretion (i.e. after accounting for 
pseudo-evolution) for haloes with masses 10 ^^ — IO^^Mq. 
However, for masses in the range 10^^ — 1 O^^M 0 , which are 
typical for LRGs, they use dark-matter only simulations 
and derive a growth of ^10% between z=0.5 to z=0 at 
scales smaller than a few hundred kpc. Together with 
pseudo-evolution, the total growth amounts to ^43%, 
co mparable to the r e sults o f iFakhouri et al.l (|2010ll . 

I Wetzel fc Nagail (|2014ll only study isolated haloes; 
LRGs cluster strongly, so the typical halo in that work 
may not be very representative for LRG s. However, no 
enviro nment selection cuts are made in iFakhouri et al.1 
(| 20 ir)ri and the derived growth rate s are similar, hence this 
seems to be unimportant. In fact, IFakhouri fc Mai (l2010ll 
measure the growth rate as a function of environment and 
find a weak trend of a decreasing growth rate towards 
denser environments, but the statistics for massive haloes 
is poor due to the low number of haloes. 

We find that the growth rates predicted from dark- 
matter only simulations are marginally consistent with 
our results. Our measurements suggest a larger growth, 
particularly towards more massive haloes. This could 
point toward the impact of baryons; particularly, AGN 
feedback may have an impact on the distribution of matter 
in massive haloes (jDuffv et al.l 120101 : IVelliscig et al.1 l2014ll 
and their accretion history. However, we need to improve 
the statistics of our measurements before we can make 
stronger claims. 

Our conclusion about the evolution of the luminosity- 
to-halo mass relation depends, however, on the validity of 
the assumption of pure passive evolution. The absence of 
a clear tilt in Figure [3] supports this view. Nonetheless, 
some residual star formation and/or mergers may also 
affect the luminosities, which could have an impact on our 
conclusions. 

The amount of star formation that is allo wed in LRGs 
is limited, based on their colour evolution (jWake et al .1 
120061 : [Maraston et al.l 1200^ . Also, spectral analyses from 
SDSS-I/H LRGs and GMASS galaxies point toward a 
very low fractio n of galaxies that ei t her form stars or hav e 
AGN activity (|Greisel et al.l l2013t I Thomas et al.1 l2013fl . 
There are indications that intermediate mass early-type 
galaxies have a low level of ongoing star f ormation (e.g: 
Kavirai et al.l 120071 : ISchawinski et al. 120071: ISalim fc Richl 
201011 . but it seems unlikely that it is sufficiently strong to 
affect our conclusions. 

Mergers also complicate a purely passive luminosity 
evolution scenario. The merging history of LRGs or mas¬ 
sive early-type galaxies has received considerable attention 


in recent vears le.e Tal et al.1 l2012l: Lonez-Saniuan et al. 

2012 

: Gabor & Dav^ 20121 iBedorf & Porteeies Zwart 

2013 

: R.uiz et al.| 20141). We focus here on the results from 

Toieiro & Percivall (l2010|l, as thev estimate the amount 


of luminosity growth in LRGs due to mergers. Their 
analysis is based on the measured luminosity function and 
clustering strength of LRGs in the SDSS in the redshift 
range 0.15 < ^ < 0.5, from which they deduce that the 
average luminosity of LRGs increases by 1.5-6% Gyr“^ due 
to mergers, depending on luminosity, such that the growth 
mainly happens for the faintest LRGs. For LRGs with 
Mrfi.i < — 22 . 8 , the evolution is consistent with passive 
evolution. In Table 4 of that work, luminosity growth rates 
from recent works are compared, suggesting that their 
val ues are fairly representati ve. 

iToieiro fc Percival (|2010ll propose two type of mergers 
that could contribute to the luminosity growth: a merger 
between an LRG and a small companion, or a merger be¬ 
tween two small companions, whose combined luminosity 
is sufficient to classify it as LRG. In the first scenario, 
the luminosities of LRGs increase towards lower redshifts; 
hence they are overestimated compared to the luminosities 
of LRGs that have evolved through pure passive evolution. 

To estimate how that may impact our derived 
luminosity-to-halo mass relations and its evolution, we 
assume that the luminosities of our three LRGs samples 
increase through mergers with 2, 4 and 6 % Gyr“^, from 
the bright to the faint bin, respectively, without affecting 
the halo masses, i.e. we assume that the mass-to-light ratio 
of the smaller companions equals zero. This provides us 
with an upper limit on the bias in the derived halo mass 
growth. We use the highest redshift slice as our reference, 
and lower the luminosities of the lower redshift slices, to 
mimic how the evolution would have looked in the absence 
of mergers. For example, between 2 = 0.6 and z = 0.2, 
approximately 3.3 Gyr passed, so we lower the luminosity 
of the Llzl bin by a factor 1.06^’^ ~ 1.2. We apply this 
factor to the average passively evolved luminosities, i.e. we 
do not apply it to the individual luminosities before the 
selection of the samples. After this adjustment, we fit the 
luminosity-to-halo mass relation again. We find that the 
retrieved slopes of the luminosity-to-halo mass relation are 
within the error bars of our nominal results. The ampli¬ 
tudes move up by ~2cr for our 0.15 < z < 0.29 redshift 
slice, and by ^Itr for our 0.29 < z < 0.43 slice. The average 
growth in halo mass becomes 116 / 135 %, consistent with our 
nominal result. We ignored here that mergers also cause 
growth in mass, hence the LRGs move diagonally rather 
than horizontally in the luminosity-to-halo mass plane. If 
the mass-to-light ratio of the smaller companions is similar 
to that of the LRG, they move along the luminosity-to-halo 
mass relation and no bias is caused; if this ratio is larger 
than that of LRGs, the actual growth in halo mass would 
be smaller than our nominal value. 

The second channel for luminosity growth through 
mergers described by Tojeiro is through the merging of two 
faint galaxies, which are individually not bright enough to 
be selected as an LRG. As the formation history differs 
from that of the typical LRG that formed at high redshift 
without much activity afterwards, the properties of their 
haloes may be different, complicating the interpretation of 
the measured trends. 

Apart from that, there are several other physical 
processes that complicate the interpretation. For exam- 
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Fig. 8. Best-fit mass-concentration relation derived for the 
0.15 < z < 0.29 LOWZ bin. The solid black line indi¬ 
cates the best-fit result and the turquoise area the 68% 
confidence intervals. The blue d o tted line shows the ref¬ 
erence relation from iDnffv et al.l (|2008ll . whose amplitude 
we fit in the halo model. The red and green dot-d ashed 
lines show the r el ations from iMaccio et al.l (|2008ll and 
iDutton fc Maccibl l|2ni4ll . respectively. The black points 
and the gr ay area are the resu l ts fro m direct fits to lensing 
data from iMandelbaum et al.l (l2008ll . 


pie, star formation and mergers may be linked in LRGs 
(jKavirai et al.l[20nll . In addition, some high-redshift LRGs 
may not become low-redshift LRGs, as they may be tidally 
disrupted or merge with other galaxies. Furthermore, 
AGN activity could trigger star formation, leading to too 
blue colours to match the LRG colour selection at lower 
redshift. The subset of LRGs whose luminosity-to-halo 
mass evolutio n is most easy to in terpret is probably the 
brightest one: iToieiro et al.f (|2012f) find that the brightest 
GMASS gal axies are most likel y to ev olve in SDSS-I/II 
LRGs, and iToieiro fc Percivriil (|20ir)ll derive that the 
brightest LRGs experience the least amount of Inminosity 
growth from mergers. 

To summarise the above: although passive evolution 
accounts for the bulk of the luminosity evolution of LRGs, 
there are several processes that complicate this picture. 
This affects our ability to select and compare LRGs and 
their progenitors at different redshifts. In order to extract 
the full information content on the evolution of LRGs that 
is contained in our measurements, one should compare with 
numerical simulations that contain all relevant physical 
processes. No luminosity evolution correction needs to be 
applied in that case, since there is no need to match the 
samples at different redshifts, as long as the selection can 
be matched in observations and simulations. The range in 
luminosity and redshift of the LRG sample studied in this 
work would form an ideal observational test case for such 
a study. 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 

Redshift 

Fig. 9. Bes t-ht normalisation o f the mass-concentration re¬ 
lation from iDuffv et al.l (l2008ll for the four redshift slices. 
Filled symbols show the results for the nominal run, open 
symbols for the miscentering run. The open symbols have 
been shifted to the right for improved visibility. 


5. Mass-concentration relation 


As discussed in Section |31 we assume a functional form 
for the mass-con centration relation, i.e. the one from 
iDuffv et al.l (l2008ll . but we allow the overall normalisation 
of this relation to vary in the fit. The resulting constraints 
for the 0.15 < z < 0.29 LO WZ bins are shown in Figure[8l 
together with the nominal IDuffv et al.l (|2008f) model. The 
constraints on /cone as a function of redshift are shown in 
Figure inland listed in Table [T] 

We hnd that the derived mass-concentration rela¬ 
tions are lower than the reference model. The largest 
discrepancy is for our lowest redshift slice, where we find 
/cone = 0.68to;o6. For the other redshi ft slices, our result s 
are about Icr below the relation from IDuffv et al.l ll2008ll . 
There is no evidence for a trend of /cone with redshift, 
although the errors of the higher redshift bins are still 
fairly large. It does suggest that the redshift scaling of this 
relation is well captured by its functional form (Equation 

113. 

In Figure 13 we also plot the relation from iMaccio et ahl 

__ __ based on dark matter only simulat ions using 

the WMAP3 cosmology (iHinshaw et ^ l2009ll . and the 
relation from iDutton fc MacciH "! 20I4I1 . who derived the 
mass-concentration relation with N-body simulations 

3 the Planck cosmology (|Planck Gollaboration et fol 
. All the relations are shown for a redshift of 0.23. 
These relations derived from dark matter only simulations 
consistently predict higher concentrations in the mass 
range that we study. 

It is possible that the low normalisation results from 
choices in our modelling. Therefore, we also derive the 
normalisation of the mass-concentration relation for the 
halo model runs where we either allow a fraction of the 
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LRGs to be miscentered, where we allow a non-zero satel¬ 
lite fraction, and for the alternative luminosity evolution 
correction runs. The results are listed in Table [5] For the 
alternative luminosity evolution run, the constraints on 
/cone do not change significantly. For the other two runs, 
the normalisations increase and become consistent with 
unity, as is shown in Figure H) By adding either satellites 
or a miscentered component, we basically shift power from 
small scales to intermediate scales. Consequently, to fit the 
same data, the leasing signal of the central LRGs needs to 
become steeper, i.e. the concentrations need to increase. 
Note that the errors on /cone increase as well, due to 
the increased freedom in the model (particularly for the 
miscentering run). 

In Figure E] we al so show the results from 
iMandelbaiim et al] (|2008ll , who derived the mass- 
concentration relation by combining lensing measurements 
for L+ type galaxies, galaxy groups traced by LRGs and 
the maxBCG cluster sample in the SDSS. The model 
used for their main results did not include a miscentering 
component, but they used a minimum scale in the ht of 
500 h“^kpc for the maxBCG sample to reduce the impact 
of miscentering. We only show the measurements that 
overlap with our range of halo masses. We find that our 
results are consistent with theirs. To wards lower masses, 
the m ass-concentration relation of iMandelbaum et ^ 
(|200i^ is in agreement with theoretical expectations, but 
at larger masses, the normalisation is lower, as can been 
seen in their Figure 5. As their lenses span a broad range 
of mass, they can fit a more flexible mass-concentration 
relation to their measurements, resulting in a broader 
confidence region compared to ours. 

The mass-co ncentration rela ti on of CMASS galaxies is 
also derived in iMivatake et al.l (|2013ll . Their results are 
very similar to ours: using their fiducial model, in which 
all galaxies are located at the centre of dark matter haloes, 
they hnd a normalisation of frn-nr. = 0.78 ± 0.12 with 
respect to the relation of iMaccio et al.l (1200711 . However, 
when they include a miscentering component identical to 
ours, the normalisation factor becomes consistent with 
unity. 

The results from IMivatake et al.l (1201,811 and the ones 
presented here suggest that the low normalisation of the 
mass-concentration relation can be explained by allowing 
for satellites and/or a miscentering component in the 
modelling, although the constraints on the miscentering 
parameters are still very weak. To understand whether 
the effect is real, one either needs to move towards larger 
data-sets whose corresponding statistical errors are much 
smaller, or alternatively, one can make use of group cata¬ 
logues with spectroscopically identihed members, to study 
the brightest group members and the satellites separately. 
The over lap between the Galax y And Mass Assembly 
(GAMA; iDriver et al.l l2009l 1201 it iRobotham et al.l 1201 if) 
and the Kilo Degree Survey (KiDS; Kuijken et ah, in prep) 
could offer an excellent dataset for such a work. 

However, it is interesting to consider alternative ex¬ 
planations, such as baryonic processes that can lower the 
concentrations. AGN feedback, for example, may affect the 
distribution of dark matter at the relevant radii, as has 
been studied in IVelliscig et al.l (I2014D . In that work, several 
baryonic processes such as cooling, supernova feedback 
and AGN feedback are implemented in hydrodynamical 
simulations to study their effect on the density profiles of 


haloes. For halo masses typical for LRGs, AGN feedback 
is the dominant process and affects the dark matter distri¬ 
bution out to several times the virial radius, such that the 
density is lower at small sca les. A similar conclusion was 
reached in iDuffv et al.l (|201Clf) : whilst cooling leads to more 
concentrated haloes, this effect is counteracted by AGN 
feedback, which could cause the concentrations to fall 
^15% below the ones from dark matter only simulations. 


6. Conclusion 

We study the evolution of the luminosity-to-halo mass 
relation of LRGs, combining SDSS photometry and 
spectroscopy with the excellent imaging data from the 
RCS2, which enables us to measure the weak lensing 
signal up to z^O.6. We use stellar population synthesis 
modelling to compute the correction to account for passive 
luminosity, which is thought to dominate the luminosity 
evolution of LRGs, enabling us to compare low-redshift 
LRGs with their predecessors at higher redshift. We split 
the LOWZ and CMASS galaxies in two redshift slices 
and three luminosity bins each, resulting in twelve lens 
samples, and fit a halo model to the lensing signals of the 
three luminosity bins simultaneously, for each redshift slice 
separately. The halo mass estimates that we obtain are 
broadly consistent with various literature results that are 
based on a variety of measurement techniques, but span 
a considerably larger combined range of luminosity and 
redshift. 

We find a typical value of ^1.4 for the slope of the 
luminosity-to-halo mass relation and no evidence that 
it changes with redshift. The amplitude of this relation, 
however, does increase significantly with redshift. We find 
that the average halo mass of LOWZ galaxies increases by 
25/^4% from (z) = 0.36 to (z) = 0.22. The halo masses 
of CMASS galaxies grows by 10 / 20 % from (z) = 0.6 to 
(z) = 0.5. If CMASS galaxies are the predecessors of the 
LOWZ ones, the total growth of LRGs from (z) = 0.6 to 
(z) = 0.22 is 80/^^%. 

This growth in halo mass is somewhat larger than what 
is expected for pure pseudo-evolution, i.e. the evolution 
in the definition of M 200 caused by the change of the 
mean background density as the universe expands, which 
by itself causes an apparent halo mass growth of ~33% 
at typical LRG masses between ^=0.6 and z=Q.2. Our 
measurements provide weak support for additional dark 
matter accretion. 

We have tested the sensitivity of these results against 
changes in the luminosity evolution and changes in the halo 
model. We find that the inferred slopes of the luminosity- 
to-halo mass relation do not significantly change. The 
amplitude of this relation decreases by at most 2cr, in the 
case where we include a satellite component in the halo 
model. The inferred average growth in halo mass of LRGs 
does not change by more than Icr. Hence for this work, 
systematic errors are likely subdominant to the statistical 
ones. For future work that uses measurements with a 
higher statistical precision, advances in the modelling of 
the luminosity evolution and of the set up of the halo 
model, are necessary to avoid biases in the results. Such 
measurements will be highly valuable to constrain the 
impact of baryonic processes on the distribution of dark 
matter, which is essential for a correct and optimized 
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exploitation of futur e cosmic shear surveys such as Euclid 
(|Laureiis et al.l[20Tlll . 

We also constrain the overall normalisation of the 
mass-concentration relation for each of the four redshift 
slices. We find that for our lowest redshift slice, the best-fit 
relation is lower than what is expected from dark matter 
only simulations. However, if we allow for miscentering or 
for the contribution from satellites in the halo model, the 
normalisation increases and becomes consistent with the 
results from dark matter-only simulations. 
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Appendix A: Luminosity evolution scatter 

The evolution corrected luminosities of galaxies become in¬ 
creasingly uncertain with redshift. This is mainly due to 
the intrinsic variation in properties of the LRGs, as illus¬ 
trated in the bottom panel of Figure[2l where the difference 
between the e-correction curves increases with redshift for 
the different models. 

To obtain a rough estimate of how intrinsic scatter in 
the {k + e)-corrected absolute magnitudes affects our re¬ 
sults, we do a simple test. We start with randomly draw¬ 
ing 10® redshifts and (fc-l-e)-corrected absolute magnitudes 
from the original distribution, i.e. the one that is shown in 
the bottom panel of Figure [31 We assume that these mag¬ 
nitudes are the intrinsic ones. Next, we assign a mass to 
each object using our nominal best-fit luminosity-to-halo 
mass relation from Table |3J Note that our conclusions do 
not sensitively depend on the choice of slope and offset. 
We compute an NFW profile for each object at 100 loga¬ 
rithmically spaced radial bins between 0.05 h^g^ Mpc and 
1 h^g^ Mpc. We stack the NFW profiles of objects that fall 
inside a lens bin as defined in Figure [31 and fit an NFW 
to the resulting profile using the mean redshift of those ob¬ 
jects. 

Next, we simulate the intrinsic scatter in absolute mag¬ 
nitude by assuming that it can be described by a Gaussian 
whose width increases as cr = 0.3z. This particular choice 
covers the range of e-correction curves of the SSP models 
that are in reasonable agreement with the colours of the 
LRGs, as shown in Figure [3] We draw a random number 
from this Gaussian and add that to the ‘intrinsic’ magni¬ 
tudes: these are the ‘observed’ magnitudes. We stack the 
NFW profiles, but now using the ‘observed’ magnitudes 
to select the lenses. Again, we fit an NFW profile to the 
stacked profile using the mean redshift of the observed ob¬ 
jects. The impact of intrinsic luminosity scatter is then esti¬ 
mated from the ratio of this ‘observed’ mass to the intrinsic 
one. We find that the mass typically changes by a few per¬ 
cent. The largest difference is found for the L3z4 bin, where 
the difference is ~10%. However, in all cases, the error is 
considerably smaller than the statistical errors. We there¬ 
fore conclude that intrinsic scatter of the luminosities can 
be safely ignored. 

Appendix B: Systematic tests 

Since this work aims at measuring the redshift dependence 
of the lensing signal, we have to be particularly careful with 
systematic effects that mimic a redshift scaling. One is¬ 
sue that requires attention is the computation of the mean 
lensing efficiency. For example, if we miss the high redshift 
tail of the source redshift distribution, the mean lensing 
efficiency would be biased increasingly low for higher lens 
redshifts, which would mimic a redshift dependence in the 
lensing signal. Motivated by the differences in mean lens¬ 
ing efficiency computed from the photometric redshift cat- 
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the different source samples, even at high redshifts, this 
suggests that the source redshift distribution is sufficiently 
accurately determined with th e latter, but not with t he for¬ 
mer (see also the discussion in iHoekstra et alll2015ll . 

Appendix C: Fitting Methodology 

We use Bayesian inference techniques to determine the pos¬ 
terior probability distribution P(A|2?) of the model param¬ 
eters A, given the data V. According to Bayes’ theorem, 

(ai, 

where P{D\X) is the likelihood of the data given the model 
parameters, P(A) is the prior probability of these parame¬ 
ters, and 

P{V) = j P{V\X)P{X)dX, (C.2) 

is called the evidence. Since we do not intend to perform 
model selection, the evidence just acts as a normalisation 
constant which needs not to be calculated. Therefore, the 
posterior distribution P{X\D) is given by 


Fig.B.l. The weighted mean of the lensing signal times 
the radius for a low-z and high-z LRG lens sample, mea¬ 
sured with different source samples as indicated in the plot. 
For both lens samples, we find consistent results, suggest¬ 
ing that our measurement pipeline is robust and does not 
introduce a redshift dependent bias. 


alogue from lllbert et al.l (120091) and lllbert et al.l (120131) at 
high lens redshifts, we perform the following test. 

We divide our source sample in six different samples, 
with magnitude cuts of 22 < nir' < 23, 23 < nir' < 24, and 
22 < nir' < 22.5, 22.5 < nir' < 23, 23 < nir' < 23.5, 
23.5 < nir' < 24. We split the LRGs in a low-redshift 
sample with 0.15 < z < 0.5 and a high-redshift one with 
0.5 < z < 0.8, and measure their lensing signals using each 
source sample. Note that for each source sample, we sepa¬ 
rately determine the shear signal ar ound random point s, the 
mean lensing efficiencies using the lllbert etld] (l2013l) pho¬ 
tometric redshift catalogues, the contamination of source 
galaxies around the lenses, and the noise bias correction. 
We measure the weighted mean of the lensing signal times 
the projected separation, as that is roughly a constant with 
radius, over the range 0.15 < r < 10 hyo^ Mpc. We show 
these values for the two LRG samples in Figure [BT] 

We find that the weighted mean of the lensing signal 
of each lens sample is consistent for the various source 
samples. Glearly, the measurements are correlated for the 
source samples that overlap in apparent magnitude, but 
the four bins of 0.5 magnitude width are more or less in¬ 
dependent (not completely as they are similarly affected 
by cosmic variance). This result strongly suggests that the 
measurement process is robust. 

W e have repeated this test using the lllbert et'ffi] 
(|2009D photometric redshift catalogue instead. For the 
0.15 < z < 0.5 LRG sample the results are consistent, but 
for the high-z one, some bins differ by 2-3tT. This points at a 
problem with the source redshift distri bution used to com - 
pute the lensing efficiencies. Since the lllbert et'ffi] (|2013D 
photometric redshift catalogue gives consistent results for 


P(A|(D) oc exp 



(G.3) 


where 


— Xesd — 


EE 

fc=i j=i 


AS(Ry-lLfc)-AS](Rj|Lfc) 
0'As(Py |Lfc) 


^ (C.4) 

AS denotes the model prediction, ctas is the correspond¬ 
ing error, Wum = 3 and Aas = 12,12,12,11 for the low 
and high redshift bins of the LOWZ and CMASS samples, 
respectively. For our fiducial model, the set of model pa¬ 
rameters is Afid = (A/inean,criogM,/conc), where: 


— Mmean is the mean halo mass for the fc-th luminosity 
bin; 

— CTiogM is the scatter of {Nc\M); 

~ /cone is the normalisation of the c{M) relation. 


When exploring model variations (see Section EH) we em- 
ploy Aniiscen — (Afid:Poffi PofF) and Agatfrac — (Affd, Af^). For 
all model parameters, we adopt a flat, sufficiently wide prior 
such that the results are not biased. 

We sample the posterior distribution of our model pa¬ 
rameters given the data using a Monte-Carlo Markov chain 
(MCMC). In particular, we implement the Metropolis- 
Hastings algorithm to construct the MCMC. At any point 
in the chain, a trial model is generated using a method 
specified below. The statistic for the trial model, Xtry: 
is calculated using Equation (IG.4I1 . This trial model is ac¬ 
cepted to be a member of the chain with a probability given 
by 


f 1-0 if xly < Xcur 

\ exp[ (Xtry Xcur)/^] if Xtry ^ Xcur 


(C.5) 


where Xcur denotes the x^ for the current model in the 
chain. We initialize the chain from a random position in 
our multi-dimensional parameter space and obtain a chain 
of ^500,000 models. We discard the first 10, 000 models 
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(the burn-in period) allowing the chain to sample from a 
more probable part of the distribution. We use this chain of 
models to estimate the confidence levels on the parameters 
and on the leasing signal, as shown in Figured 

A proper choice of the proposal distribution is very im¬ 
portant in order to achieve fast convergence and a reason¬ 
able acceptance rate for the trial models. The posterior dis¬ 
tribution in a multi-dimensional parameter space, such as 
the one we are dealing with, will have degeneracies and 
in general can be very difficult to sample from. We have 
adopted the following strategy to overcome these difficul¬ 
ties. During the first half of the burn-in stage, we chose 
an independent Gaussian proposal distribution for every 
model parameter, as is common for the Metropolis-Hastings 
algorithm. Half-way through the burn-in stage, we perform 
a Fisher information matrix analysis at the best-fit model 
found thus far. The Fisher information matrix, given by 


F- ■ — 


\nC 

dXidXj 


(C.6) 


is a Np X Np symmetric matrix, where Np denotes the num¬ 
ber of parameters in our model, and C oc e~^ is the likeli¬ 
hood. The inverse of the Fisher matrix gives the covariance 
matrix, C, of the posterior constraints on the model param¬ 
eters. More importantly, the eigenvectors of the covariance 
matrix are an excellent guide to the degeneracies in the pos¬ 
terior distribution, and the corresponding eigenvalues set a 
scale for how wide the posterior ought to be in a given di¬ 
rection. Therefore, for the second half of the burn-in period, 
we utilize this information and use a proposal distribution 
which is a multi-variate Gaussian centered at the current 
value of the parameters and with a covariance given by the 
aforementioned matrix. In practice, the trial model (Atry) 
can be generated from the current model (Acur) using 


A 


try 


— Acur T C Ax, 


(C.7) 


where x is a vector consisting of Np standard normal de¬ 
viates, the matrix A is such that AA^ = C, and ^ is a 
parameter that we have chosen to achieve an average ac¬ 
ceptance rate of ^30%. We repeat the Fisher matrix anal¬ 
ysis once again at the end of the burn-in period (using the 
best-fit model found thus far) and use the covariance ma¬ 
trix to define our proposal distribution to be used for the 
MCMC. This strategy has proven to be extremely efficient 
in s ampling posterior di s tribut ions for similar studies (see 
e.g. ICacciato et aDl2013L[2014f) . 


Table D.l. Properties of the lens bins (after fc-correction) 



Mr 

Nlens 

{z) 

(Lr) 

Meft 

fconc 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

0.15<2<0.29 (LOWZ) 



2.60l°:®^ 


Llzl 

[-21.8,-21.2] 

1110 

0.204 

0.68 


L2zl 

[-22.4,-21.8] 

3784 

0.225 

1.04 

4.48l°;®« 

0 . 781^9 

L3zl 

[-22.8,-22.4] 

954 

0.234 

1.59 

n yo + 1-17 

o. i J-i.os 


0.29<2<0.43 (LOWZ) 



o 1 '7+0.50 
O.i-i -0.46 


L2z2 

[-22.4,-21.8] 

4725 

0.351 

1.10 


L3z2 

[-22.8,-22.4] 

2901 

0.368 

1.62 

c 70+0.83 

0. < 0_Q 7g 

0.85l°:1^2 

L4z2 

[-23.2,-22.8] 

981 

0.374 

2.31 

8.64l];f5 


0.43<2<0.55 (CMASS) 



2.101°+ 


L2z3 

[-22.4,-21.8] 

7798 

0.496 

1.09 


L3z3 

[-22.8,-22.4] 

4197 

0.506 

1.62 

2.481°;°^ 

0.85l°;l? 

L4z3 

[-23.2,-22.8] 

1646 

0.506 

2.31 

7.18ll:°t 


0.55<z<0.70 (CMASS) 



1.08l°;^« 


L2z4 

[-22.4,-21.8] 

2532 

0.584 

1.17 


L3z4 

[-22.8,-22.4] 

4357 

0.599 

1.67 

9 7Q+0.89 
< ^-0.77 

0.69l°;?^ 

L4z4 

[-23.2,-22.8] 

3365 

0.615 

2.35 

3.97+;i® 


L5z4 

[-23.6,-23.2] 

1314 

0.625 

3.37 

6.98+;?4 



Notes. (1) absolute magnitude range (after fc-correction); 
(2) number of lenses; (3) mean redshift; (4) mean lu¬ 
minosity [10^^Lq] (after fc-correction); (5) best-fit halo 
mass (6) best-fit normalisation of the mass- 

concentration relation. 


be more straightforward. 

Hence we divide the lens sample into bins of fc-corrected 
absolute magnitude and redshift. Details of the lens samples 
can be found in Table iDTl The selection is also illustrated 
in FigureO For each subsample, we stack the leasing signal 
of all the lenses inside that bin and show that in Figure IDTI 
We fit the halo model using the same set-up as before and 
show the best-fit models, together with the model uncer¬ 
tainties, in Figure lD.il The corresponding effective masses 
and normalisations of the mass-concentration relation can 
be found in Table IDTI 


Appendix D: k-corrected only results 

In this work, we have attempted to account for the evolu¬ 
tion of the luminosities of LRGs by applying a luminosity 
evolution correction. Unfortunately, this correction is un¬ 
certain and will also add scatter to the corrected luminosi¬ 
ties as LRGs cover a range of intrinsic properties such as 
formation age and metallicity, which complicates the in¬ 
terpretation of the results. Therefore, it is also interesting 
to apply only the fc-correction to the absolute magnitudes. 
The resulting luminosities are closer to the real luminosi¬ 
ties, i.e. they have smaller scatter. Interpreting any trend 
in the luminosity-to-halo mass relation will be harder as 
both the luminosity and the halo mass may evolve simulta¬ 
neously, but comparison to for example simulations should 
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1-22.4,-21.8] 



0.15<z<0.29 




Projected separation [/lyo^Mpc] 



Fig. D.l. The lensing signal AS of LOWZ (top two rows) and CMASS (bottom two rows) lenses as a function of 
projected separation for the four luminosity bins (after the fc-correction is applied). The solid red lines show the best-fit 
halo model, the orange and yellow regions the 1 and 2a model uncertainty, respectively. We fit the signal on scales 
between 0.05 and 2 Mpc. 
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